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An analytic column density profile to fit prestellar cores 



o 
o 

(N 



o 

(N 
> 

l> 

O 
On 
O 



Wolf B. Dapp^* and Shantanu Basu^, 

^Department of Physics and Astronomy, The University of Western Ontario, London, Ontario, N6A 3K7, Canada 



ABSTRACT 

We present a new analytical three-parameter formula to fit observed column den- 
sity profiles of prestellar cores, ft represents a line-of-sight integral through a spher- 
ically symmetric or disc-like isothermal cloud. The underlying model resembles the 
Bonnor-Ebert model in that it features a flat central region leading into a power-law 
decline oc in density, and a well-defined outer radius. However, we do not assume 
that the cloud is in equilibrium, and can instead make qualitative statements about 
its dynamical state (expansion, equilibrium, collapse) using the size of the flat region 
as a proxy. Instead of having temperature as a fitting parameter, our model includes 
it as input, and thus avoids possible inconsistencies. It is significantly easier to fit to 
observational data than the Bonnor-Ebert sphere. We apply this model to L1689B 
and B68. We show that L1689B cannot be in equilibrium but instead appears to be 
collapsing, while our model verifies that B68 is not far from being a hydrostatic object. 

Key Vifords: ISM: clouds - ISM: globules - stars: formation 



1 INTRODUCTION 

' In recent years, observational advances have made it 
] possible to measure column density profiles in prestel- 
. lar co res. Various methods are beins u s ed in the liter- 
' ature. IWard- Thompson. Motte. fc Andrj (Il999l ) meas ured 
, mm continuum emission, whil e iBacmann et aL f ll2000l) uti- 
lized mid-IR absorption, and lAlves. Lada. fc Ladal (|200ll ) 
measured dust extinction and reddening of the light of back- 
,, ground stars in the near-IR. A fourth method is to use flu x 
' measurements in optically thin lines l|Tafalla et al.l |2002| ) . 
Most of the column density profiles measured in these ways 
show certain prominent common features: a central fiat re- 
gion, followed by a power-law decline. Some cores addition- 
ally exhibit signatures of steepening of the profile, while 
some show a more-or-less smooth merger to some back - 
ground value of the column density l|Bacmann et al.ll2000l ). 

One ni o del o f ten used to fit to such cor es (e.g., 
lEvans et all I2001I : iTeixeira. Lada. fc Alvej l2005h is the 
Bonnor-Ebert sphere (hereafter BE sphere). This model as- 
sumes an isothermal sphere in equilibrium, acted upon only 
by gravity and thermal pressure, a nd bo und by some exter- 
nal pressure (Bonnor 1956; Ebcr Hl955l ). In the central re- 
gion its density profile is fiat with density « = (r = 0). 
The size of this region is proportional to the Jeans length, 
Rj oc Cs/\/GQc- It transitions into a power- law decline 
that approaches the Singular Isothermal Sphere (SIS) with 
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esis = cl (27rGr2) ^ for large radu (|Shulll977l ). The cloud is 
finally cut off at some finite radius, where external pressure 
forces match the internal forces. The BE model invokes the 
dimensionless radius parameter ^ = r^^A-KGQc/ca- The value 
^crit = 6.45 marks a dividing line such that clouds trun- 
cated at dimensionless radius > Ccrit are in an unstable 
equilibrium, and clouds with < ^crit are stable. 

The column density profile of the BE sphere, obtained 
by integrating the volume density numerically al ong a line of 
sight , matches well with some observations (B68: lAlves et all 
120011 ). Depending on the parameters, the power-law decline 
oc r~^ in column density can be more or less pronounced, or 
even almost completely absent. The profile steepens at the 
edge because the line of sight through the sphere becomes 
shorter. This is a geometric effect which is present in all 
truncated models, unless the density increases sufficiently 
fast with radius. 

While physically motivated and reproducing features of 
several observed column density profiles, the BE model has 
shortcomings. The most important is the key assumption 
of e quilibrium. In fact, m ost fits are found to be supercrit- 
ical (|Teixeira et al.l [20051 ). representing unstable equilibria. 
These states are not expected to exist in reality, as any per- 
turbation will send them to immediate collapse. On a more 
practical side, the procedure of fitting the BE model to ob- 
servations is quite involved. The volume density is only avail- 
able as a numerical solution, which then needs to be inte- 
grated (again numerically) to calculate the column density. 
Sometimes the fit demands temperatures well above those 
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meas u red for the centres o f prestellar cores llBacmann et al 



200d: lAndre et al] l2003l : iKirk. Ward-Thompson, fc Andre 



20051 ). and predicted by de tailed models of the thermo- 



dyna mics within the core (iGalli. Walmslev. fc Goncalved 
12002! ') . Finally, most of the observed cores are deeply embed- 
ded within their parent clouds, and the source of a bounding 
exter nal pressure is no t obvious. The often-cited example 
B68 (|Alves et al.ll200ll ) is an exception, as it is thought to 
lie within a hot HII region. 

We stress that the generic features of the BE den- 
sity profile are not unique to equilibrium situations (see 
iKandori et al] |2005| ). The flat region with adjacent r 
density profile appears also in solutions of the hydro- 
dynamic equ ations for gravitationally collapsing objects 
llLarsorJ Il969l). as long as pressure is not completely neg- 
ligible (|Shu. Adams, fc Lizanolll987l ). The pressure gradient 
then establishes a region where it nearly balances gravity. 
Here, the density is nearly constant on the scale of the lo- 
cal Jeans length, shrinking in size over time as the density 
increases. Inverted, this requires that g cc in the outer 
profi le that is left behind outside the central region fsee lBasul 
1 19971 ). 

There is no reason why the above two features (flat re- 
gion with size of Jeans length and adjacent density pro- 
file) sho uld be present in a non-self-gravitating cloud. How- 
ever, iBallesteros-Paredes. Klessen. fc Vazguez-Semadenil 
l|2003t) find that convergent turbulent fiows, with and 
even without self-gravity (which are expected to have 
very different volume density profiles), nevertheless yield 
column density profiles that resemble those of a BE sphere. 
There are three reasons for this seemingly surprising 
result. The first one is the effect of smoothing the data, 
by angle- averaging, and also by integration along the line 
of sight. Secondly, the BE sphere leaves the modeller the 
freedom to fit the size of the central flat region by v arying 
the temperature. Indeed, iBallesteros-Paredes et aU (|2003 ) 
flt their simulated cores with BE temperatures in the range 
of 5 — 60 K, despite their models being strictly isothermal 
at T = 11.3 K. Finally, the position of the outer radius 
cutoff is somewhat arbitrary. This often leaves large parts 
of profiles unfitted. In this paper, we argue that fitting a 
prestellar core profile at a set temperature does still allow 
one to distinguish between different models of internal 
structure. 

We propose an analytic density profile reprodu cing the 
characteristics of not only isothermal equilibria (|Bonnorl 
1 19561 : lEberdll955l'). b ut als o non-equilibrium collapse so- 



lutions (e.g., LarsonI Il969l ) , and many observed profiles 



(|Bacmann et alj 2000). Within the margins of uncertainty 
it fits the observations as well as the BE model does. How- 
ever, it possesses a closed-form expression for the column 
density, and is therefore very easy to fit. Furthermore, the 
temperature can be an input parameter instead of a fitting 
parameter, so that the model avoids some of the possible in- 
consistencies of fitting the BE model to either observations 
or simulation results. We can use our model to make some 
inference about the dynamical state of the core. 

This paper is organized the following way: Section [2] 
describes the spherical model and its parameters, while Sec- 
tion [3] presents a corresponding model for intrinsically flat- 
tened objects. In Section [l] we apply our model to B68 and 
L1689B, and we summarize our results in Section [5] 




Figure 1. Schematic illustration of a cut through a spherical 
cloud of radius R. The observer is positioned along the direction 
of the coordinate s, and measures an integrated column density 
S as a function of the offset x. 



2 SPHERICAL GEOMETRY 
2.1 Basic model 

The characteristics found in observed column density pro- 
flies and theoretically both in equilibrium and collapse solu- 
tions can be parametrized by a volume density more generic 
than the BE sphere. We propose using the profile 



g{r) = 







r>R, 



(1) 



characterized by a central volume density Qc and truncated 
at some radius R. The parameter a fits the size of the flat 
region in terms of the Jeans length, and is given by 



a = k- 



(2) 



where G is the gravitational constant and k is a cons t ant o f 
proportio nality. This pro flle is also mentioned in iKind (|l962l ) 
and Tafall a et alj l|2002f ). The temperature T, which can be 
constrained observationally, is not used as a fltting param- 
eter. It enters through the value of the isothermal sound 
speed Cs — ^/ksT/Jt. The Boltzmann constant is denoted 
by kB, and fi = 2.33 mn is the mean mass of a particle, 
where ttih is the mass of a hydrogen atom. We assume a 
10% number fraction of helium. 

The column density is derived by integrating the volume 
density along a line of sight through the sphere: 



E{x) = 2 



Q (s) ds 



g (r) rdr 



(3) 



and 



where we have used the transformation s — \/r^ - 
hence ds = rdr/y/r'^ — x'^. Fig. [T]deflnes the quantities ap- 
pearing in this derivation. 

Inserting equation ([1} into equation Q, we flnd that the 
integral is analytically tractable. The closed-form expression 
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Figure 2. Plots of equation JSj for varying c = R/a. All cases 
have identical values for a. The central column densities are cho- 
sen to be different in order for the plots not to overlap. The solid 
line illustrates the three major characteristics of this function: 
the fiat region in the centre of size a, the adjacent power-law 
decline oc x~^, a steepening at the edge for geometric reasons 
(see text), out to the cut-off radius R. In this case, c = 50. The 
dashed-line shows a model dominated by the power-law regime, 
with c > 1000. R is chosen so that the influence of the cut-off 
is negligible. The long-dashed line finally demonstrates that the 
power-law can be suppressed entirely if R > a. In this case, c = 5. 



for the column density is then 



E(x) = 



: arctan 



(4) 



This can be re- written in terms of the central column density 
by introducing the ratio c = R/a, called the concentration 
in iKind (|l96a ). With E = 0) = Ec = 2agc arctan (c) we 
find 



arctan 



{x/aY 



1 + {x/aY 



arctan (c) 



(5) 



Fig. [2] demonstrates that this profile possesses generic 
features that it shares with observations, collapse solutions, 
and with the BE model: a fiat central region, a power-law 
decline, and steepening at the edge. We note that the effect 
of the boundary is exclusively contained in the factor in 
square brackets. 

The quantity c determines the size of the region de- 
scribed by the power-law. If it is large, there is a pronounced 
power-law, whereas if c approaches unity, the cut-off already 
dominates near the fiat region and inhibits the power-law. 

The three parameters to fit are (i) the outer radius R, 
(ii) the central column density Ec, and (iii) the size of the 
fiat region a. The latter contains the product of k and VT, as 
shown in equation ([2]). If the temperature is pre-determined, 
a only depends on the value of k. This then allows for a sta- 
bility assessment, as demonstrated in the following section. 



2.2 Model parameters 

The pressure in an isothermal system is given by P = c^g, 
and hence the pressure gradient for the spherically symmet- 
ric volume density profile of equation is 



dP 
dr 



idg 
' dr 



a2 [l + ir/af]^' 



(6) 



In the inner regions, where r <^ a, 



(7) 



dP _ -2clQcr 
dr 

This demonstrates a simple point: the larger the fiat re- 
gion a, the smaller the pressure gradient in that region. The 
minimal value for a is reached in an equilibrium situation. 
A larger a would result in a smaller pressure gradient, grav- 
ity would win out, and collapse would ensue. Conversely, 
for a smaller than its equilibrium value the pressure force 
would dominate, causing expansion. However, that case is 
not probable since a profile with g Qcr~^ mandates a strong 
gravitational infiuence. 

Another well-known model besides the SIS is the so- 
called Larson-Penston solution (LP solution), a self-similar 
spherical collapse soluti on, assuming a homogeneo us ini- 
tial density distribution l|Larsonl [19691 : |Penstonlll969l ). This 
highly dynamical model does not assume equilibrium, and 
asymptotically reaches a density profile for which 



f>LP = 4.4 gsis 



(8) 



everywhere. We now want to find an expression for the pa- 
rameter a in the SIS and the LP models. In the outer regions, 
equation ([T} becomes 



£>(r » a) 



gca 



(9) 



■2 2nGr^ ' 

where k = 2Tvk^ — 1 for the SIS (equilibrium) and k = 4.4 
for the asymptotic LP collapse solution. Therefore 



0.399 for SIS, 
0.837 for LP. 



(10) 



This has a direct application to observed prestellar cores. If 
an estimate for the core temperature is available, equation 
(O can be fit to the column density profile, with fc as a 
fitting parameter. If A; ~ 0.4, the system can be considered 
to be in equilibrium. If k is significantly larger than 0.4, and 
closer to 1, the cloud under scrutiny is collapsing. 

Fig. |3] shows our model with a chosen to converge to 
either the SIS or the asymptotic LP solution. We also plot 
the equilibrium BE model. In order to avoid boundary ef- 



fects, the outer radius was moved to 



The 



BE profile does not join onto the SIS profile right away. 
It overshoots, then steepens and asymptotically approaches 
the SIS profile at large radii. 

We may ask which value of k is most appropriate to 
model a critical BE sphere. The overshoot effectively pro- 
duces a larger fiat region if the solution is truncated at 
= Ccrit- The resultant a has to be larger than for the 
SIS. Having fixed the temperature this can only be achieved 
by increasing k, as we show in Fig. ID Profile (a) shows a 
subcritical BE sphere, i.e. a stable equilibrium solution. In 
this case, there is no discernible power-law region and the 
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Figure 3. Comparison between the BE sphere (dashed line) 
and our model (solid lines) for identical central temperature and 
cut-ofT radius. The top line is our model with k = 0.837, for 
which the profile approaches the asymptotic LP solution. The 
lower line shows k = 0.399, for which the SIS profile is reached. 
The dashed line represents the BE model, which also asymptoti- 
cally approaches the SIS. 



central flat region makes up a significant portion of the to- 
tal radius. The cut-off becomes dominant just outside the 
flat region, and the best-flt dimensionless dynamics param- 
eter is fc = 0.63. Note that the density contrast is less 
than an e-folding, which means a mean column density en- 
hancement of < 40% over the background. Such an object 
woul d not be observationally charac terized as a prestellar 
core l|Andre. Basu. fc I nutsuka"2009Y 

Profile (c) in Fig.|4]shows the opposite situation, where 
the power-law portion of the profile is fit, and k = 0.46 is 
closer to the equilibrium value. As shown in Fig. |3l the BE 
model initially has a power-law index steeper than —2. This 
is the reason why the profiles do not match up as well as for 
the other cases. Profile (b) finally shows a critical BE sphere. 
A density contrast greater than Qc/ Qs = 14.0 would charac- 
terize an unstable equilibrium, initiating collapse upon the 
smallest of perturbations. The two profiles are very simi- 
lar, only minor deviations are discernible. In order to fit our 
model to the overshoot of the BE profile over the SIS, a has 
to be larger, which can only be achieved by fc = 0.54. We 
demonstrate this effect on the examples of B68 and L1689B 
in Sections 14.11 and 14.21 as a consistency check, we insert 
the derived BE temperatures into our model, and indeed 
find fc = 0.57 and k = 0.56 respectively. Note that the BE 
model is a valid model in absence of any effects other than 
gravity and thermal pressure, and thus we conclude that in 
this case, not the asymptotic value of fc = 0.4 but rather 
the critical BE value of fc = 0.54 is relevant for the stability 
assessment. We stress that if k is significantly larger than 
that, it is strongly indicative of collapse. 

Since the BE model effectively has fc ~ 0.5 and coUaps- 
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Figure 4. Fit of the integrated column density of various BE 
spheres (dashed lines) with our model (solid lines), (a) is a sub- 
critical (stable) model, while (b) shows a critical BE sphere. In 
that case, k = 0.54 (see text for explanation), (c) presents a fit 
to a highly supercritical (unstable) BE model. 



ing clouds have fc ~ 1.0, trying to force a BE model fit to 
such a cloud r esults in spuriou s high er temperatures. This 
can be seen in iBacmann et all Hooa). A factor of 2 in the 
dimensionless dynamics parameter fc needs to be compen- 
sated by a factor of 4 in temperature when fitting a BE 
equilibrium state. 



2.3 Total mass 

Both our model and the BE sphere have a well-defined outer 
radius, and the total mass depends on it. For the profile of 
equation ([T}, the total mass has a closed-form expression. 
Evaluating 4Trr^g{r) dr yields 



M = 47r£ica 



arctan 



R 



(11) 



There is evidence for some power-law beha viour in most ob- 
served profiles (e.g., iBacmann et al.ll2000l ). In such a situ- 
ation, the value of R/a must be at least > 10 (cf. Fig. [2]). 
Then, the contribution of arctan (c) is quite close to its lim- 
iting value of 7r/2, and the total mass is given by 



M ^ Angca^R = A-Kk^RcJG 



(12) 



with < 15% error. The equivalent expression for the BE 
model is (derived assuming an SIS profile at large radii) 



A/be ~ 2RcJG. 



(13) 



Both expressions for the masses are linearly dependent 
on the unknown radius R. For a given (measured) temper- 
ature, our model additionally invokes the dimensionless dy- 
namics parameter fc, which is well constrained by the size of 
the flat region. The BE mass instead depends on the tem- 
perature of the fit, which we show in Sections 12. 21 and 14. 21 to 
potentially be a very poor value, and prone to overestima- 
tion by as much as a factor of 4. 
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3 DISC GEOMETRY 
3.1 Basic model 

Most observed cores do not appear circular in projection. A 
corresponding analysis as in Section 12.21 can be done for a 
disc-like geometry as well. Here, a generic face-on column 
density profile is 



ac ^ ^Jl + {r/af r sC i?, 



(14) 







r>R. 



This profile provides a goo d fit to the column density of 
collapsing fiattened clouds (jBasul Il997h , even though the 
models in that paper are not truncated. 

Assuming vertical hydrostatic equilibrium and ignoring 
the effect of external pressure, the volume density is propor- 
tional to the square of the column density 



2 



and the density accordingly is given by 



ttG 



2ci 1 -I- (rjaf 



(15) 



(16) 



The assumption of vertical hydrostatic equilibrium is well 
justified by simulations when some source of flattening 
is pr esent, such as rot ation jNarita. Havashi. fc Mivamal 
1 19841 ). or magnetic fields (jPiedler fc MouschoviaJ 19931 ). Ab- 
breviating -KGa"^/ (2Cs) = Qc, the volume density profile has 
exactly the same form as for the spherical case (with a dif- 
ferent constant). 



3.2 Model parameters 

Integrating q (r) through the disc viewed edge-on, the col- 
umn density versus the ofi^set x has the same functional form 
as derived above for the spherical case in equation ([5}. The 
only difference is that the central column density is now 
given by Ec = tt/cctc arctan (c), since for a thin disc, the re- 
lation between a and the Jeans length is 



k 



(17) 



This similarity means that both a flattened and a spher- 
ical object (see Section [2| can be flt with the same formula. 
Note that o (r) is the face-on column density, whereas E [x) 
is its edge-on counterpart. 

We can compare the (face-on) generic column density 
profile to one of an equilibrium solution, in order to get the 
minimal value for a, as done for the spherical case above. 
Here, the appropriate profile is the singular isothermal disc, 
characterized by 



- 2^- 

Following the same steps as in Section [2.21 the large-radius 
asymptote of equation (|14p is 



cr (r > a) 



(19) 



yielding a = /(2TvGac). Comparing this expression to 
equ ation l|17p. we find k = l /27r = 0.160 for equilibrium. 
ISaigo fc Hanawal (| 19981 ) show that for an isothermal 



disc-like cloud during runaway collapse, the self-similar col- 
umn density profile is given by 3.61 asm, analogous to how 
the LP solution is over-dense compared to the SIS. This 
provides a good estimate for the value of the dimensionless 
dynamics parameter for a dynamically collapsing disc-like 
cloud, k = 0.57. Flattened cloud cores best fit with values 
for the dimensionless dynamics parameter k m 0.2 are there- 
fore close to equilibrium, while k ~ 0.6 is strongly indicative 
of dynamical collapse. 



3.2.1 Effect of magnetic fields 

Inclusion of magnetic fields into a spherical model poses a 
problem: the magnetic field lines cannot be arranged in a 
spherically symmetric way. There has to be a preferred di- 
rection. However, in disc-like geometry, the axis perpendicu- 
lar to the disc is preferred already, and can be chosen as the 
axis of orientation of a large-scale magnetic field. Starting 
with initially straight field lines through a collapsing astro- 
physical gas sphere, flattening ensues along this preferred 
direction. Matter can contract significantly along the field 
lines, but the Lo rentz force will imp ede motions perpendic- 
ular to the field (|Mouschoviaslll976l ). 

In the limit of a magnetic field much st ronge r 
the ambient cloud, iBasul l)l997h . 



than the field of t he a mbie nt cloud, iBasi] 
iNakamura fc Hanawal l| 19971 ') and IShu fc Lil l| 19971 ') showed 
that the contribution of the magnetic field can be folded 
into a force calculation of the collapse of a magnetized disc 
by magnetic pressure modifying the effective sound speed, 
and magnetic tension changing the effectiv e gravitational 
constant. Then we can write fsee lBasulll998l ) 



a = k 



cl (1 + 2^l-^) 



(20) 



We use the standard defi nition of the mass-to-flux ra tio in 
units of the critical value l|Nakano fc Nakamura|[l973 ): 



B. (r) 



2ttVG. 



(21) 



Collapse requires /i > 1, otherwise the magnetic forces will 
dominate over gravity and stabilize the cloud, and evolution 
can only happen on long time scales by ambipolar diffusion 
(see, e.g.. lShu et al]|l987l ). 

We can absorb the modification into an effective k, 
writing it as fcetr = A; (l -I- 2fj,~^^ (l — • The value 

k — 0.160 as derived for the non-magnetic case remains the 
minimal value, while feeff can exceed it by a factor of 1.1 to 
2 for /i « 2 — 5, which is a reasonable estimate for supercriti- 
cal co res (jCiolek fc Mouschovias|[l994l : lBasu fc Mouschoviad 
11994!). This means k^s < 0.3, which leaves it still a factor 
of 2 smaller than the value for the dimensionless dynamics 
parameter for dynamical collapse (fc ~ 0.6, as shown above), 
and thus clearly distinguishable. 



3.2.2 Effect of rotation 

Similar to the problem for magnetic fields, rotation cannot 
be considered in strict spherical symmetry. In fact, rotation 
can be the cause of disc-like geometry, as a collapsing rotat- 
ing spherical cloud will ffatten along the rotation axis before 
contracting in the radial direction. 
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There exist self-similar solutions for the collapse of a ro - 
tating thin disc l|Narita et al. II 1984 ISaigo fc Hanawalll998h . 
Following the discussion in " iBasul l|l99'it . one can express 
the additional effect of rotation as an effective acceleration 
acft = aT (1 + ac/ax), where st is the thermal acceleration, 
while ac denotes the centrifugal acceleration. Assuming the 
column density profile of equation (|14p and proportionality 
between the specific angular momentum and enclosed mass, 
this leads to 



ac/flT ~ 3 X 10 



(22) 



This number is computed for a background rotation rate of 
the molecular cloud Q,c = 10~^* rad s~^, a central column 
number density TVc — 10^^ cm~^, and a temperature T = 
10 K. Hence, we find that the effective radial acceleration 
opposing gravity exceeds the thermal acceleration by less 
than 1%. This shows that, unlike magnetic forces, rotation 
does not significantly modify the size of the flat region and 
hence the dimensionless dynamics parameter k. 



4 APPLICATIONS 

We fit both the BE model and our model to observational 
data by determining the best-fitting parameters using a 
standard Levenberg-Marquardt least-squares minimization 
algorithm based on MINPACK. 



4.1 B68 

We fit our model to B68, the prime example in the liter- 
ature for an extraordinarily good fit of the BE model to a 
(angular-averaged) column density p rofile measured in near- 
IR dust extinction ijAlves elaLlliool). B68 is an isolated Bok 
globul e, a small dark clo ud which has been studied exten- 
sively. lAlves et al.1 (|200lh assume a distance of 125 pc and 
quote a BE mass of 2.1 M©, a temperature of T = 16 K, 
an outer radius of 12, 500 au, and a dimensio i iless outer 
radius ^ = 6.9 ± 0.2. iHotzel. Hariu. fc Juvelal l|2002h and 
iLai et all (|2003[ l more recently updated some of these val- 
ues by measuring the temperature to be T = 10 ± 1.2 K 
and T = 11 ± 0.9 K, respectively, and estimated the cloud 
to be closer by about 25 pc, placing it at the near side of 
the Ophiuchus complex. This reduces the outer radius and 
decreases the BE mass to « 1 M©, but does not change the 
value for ^a, as that is determined by the shape of the profile . 

We perform a similar analysis as lAlves et al.1 (|200lf ). 
We calculate a BE sphere, vary the temperature, the di- 
mensionless and physical outer radii and R, and fit the 
line-of-sight integral to the observational data, assuming a 
distance of 100 pc. This procedure yields a BE mass of 
M = 1.17 Mq, a BE model temperature of T = 11.1 K, and 
a central number density of ric — 2.3 x 10^ cm~^. The outer 
radius of the best BE fit is 10, 680 au, and the dimensionless 
outer radius = 7.0. The best fit of our model to the same 
data yields a total mass of M — 1.2 A/0, R = 10,420 au, 
and Tie — 2.7 x 10^ cm~^. The size of the fiat region is 
a = 2,830 au. Assuming a temperature of T = 11 K, as 
in the observations mentioned above, this corresponds to 
k = 0.57. We showed in Section r2.2l that a critical BE sphere 
is fit with k — 0.54, and we conclude that the internal struc- 
ture of B68 may indeed very closely resemble a critical BE 
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Figure 5. Best fit of the BE model (da shed line) and our model 
(solid line) to the column density of B68 l lAlves et al.ll20o"lh . Both 
fits follow the data points extremely closely and yield very similar 
results for the total mass (see text). Like the BE model, our model 
fit is consistent with an equilibrium cloud. 



sphere. Fig. [5] shows that the best fits to B68 of the BE 
sphere and our model differ by very little over the whole 
range of data. 

The BE analysis is much more involved and computa- 
tionally expensive than our model. It necessitates a numeri- 
cal solution of the Lane-Emden equation. This ordinary dif- 
ferential equation (ODE) underlies the BE model, and does 
not have a general analytical solution. The numerical so- 
lution is truncated at the dimensionless radius and con- 
verted into physical units. Only then can one integrate along 
lines of sight through the solution numerically, yielding the 
column density to be compared with observational data. In 
contrast, fitting our model to a dataset requires less than a 
dozen lines of code. 



4.2 L1689B 

Another prestellar core to which we apply our model is 
L1689B. Its central column density is much greater and it is 
more extended than B68. Furthermore, it is not isolated, but 
rather embedded within the larger complex o f the Ophiuchus 
molec ular cloud. Its profile was measured bv lBacmann et ahl 
using m i d-IR absorption observations and updated 
bv lAndre et al.l (|2003l ). We use the East- West profile, which 
shows evidence of steepening in the outer regions. 

The BE fit to the data in Fig.[6]yields = 11.2 (super- 
critical, unstaW«_£ore)_aiid a mass of M = 9.3 Mq, similar 
to what lAndre et all (|2003l ') found. Fitting our model results 
in M = 9.1 Mq. However, in contrast to the case of B68, 
the temperatures of the two fits are vastly different. The 
BE model requires a temperature T = 40 K, whereas our 
model can achieve the same quality of fit with a temperature 
T = 10 K. We fit the size of the fiat region a = 3, 600 au 
which is a combination of the temperature and k. A temper- 
ature of T = 10 K mandates k — 1.1, which is very far from 
the value required for equilibrium, and even surpasses the 
value k = 0.837 achieved in the highly dynamical asymptotic 
LP solution. As discussed in Section [2.21 this is a strong in- 
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Figure 6. Best fit of the BE model (dashed Hne) and our model 
(solid line) t o the column density of L1689B, as measured by 
lAndre et al. I l l2003t) . A gain, both fits are very similar and yield 
observationally almost indistinguishable results. However, the BE 
model requires a relatively high temperature of T = 40 K, while 
our model can fit the data with T = 10 K. The dimensionless 
dynamics parameter k = 1.1, which means that the cloud is col- 
lapsing. 

dication that L1689B cannot be in equilibrium, but must 
instead be collapsing. 

Again, our model fits this dynamically evolving object 
just as well as it fit the presumably almost static B68. It 
can do so without introducing inconsistent temperatures, 
and with much less computational and coding eflort. 



5 SUMMARY AND CONCLUSIONS 

We have introduced an analytic profile for the integrated 
line-of-sight column density of an isothermal spherical or 
flattened cloud. This cloud can either be in equilibrium or 
in a state of dynamical collapse. Our model is very simple 
to calculate compared to the BE model, and a few lines of 
code suffice to find a best-fit set of parameters. Another ad- 
vantage lies in its ability to also encompass non-equilibrium 
states. The dimensionless dynamics parameter k allows one 
to assess whether a cloud is near equilibrium or vigorously 
collapsing. At the same time, our model does not produce in- 
consistencies like the BE model regarding the object's tem- 
perature (which can be treated as a constrained quantity 
and not a free parameter). 

Our model can be applied both to spherical and to flat- 
tened clouds, and in both cases yields the same functional 
form for the column density E (x). It fits the size of the cen- 
tral fiat region and allows the modeller to adjust either k or 
the temperature to match it. 

Results of fitting our model to B68 show that it is indeed 
a near-equilibrium cloud, with parameters very similar to 
the best-fit BE model. For L1689B, our model avoids the 
need for a high temperature (the BE model requires T = 
40 K) since it can be interpreted as having a temperature 
of T = 10 K but being in a state of dynamical collapse. 
This finding is confirmed by the detection of infall motions 
for L1689B by means of the shape of the line profiles in 
optically-thick molecular transitions (jBacmann et al.ll2000l : 



iLee. Myers, fc Plumdiiooi ). In the future, our model can be 
applied to many other prestellar cores. 
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